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wave pulse response is generated twice (h 3 (k) and ^(k)), 
and square-wave pulse responses (h 3 (k) and h 4 (k)) are 
folded together so that the elements 9, 10, 5a, 11, 12 
generate the triangular pulse response hj(k) shown in FIG. 
5 on the whole. By sending the output signal of the second 5 
differentiator 12 to the averaging filter 13, there is a folding 
with the Dirac pulses of the pulse response of the averaging 
filter 13, so that the triangular pulse response hj(k) is added 
to the triangular sequence offset by the period Vz T fl2 (offset 
by 40 sampling times in the embodiment illustrated in FIGS. 10 
5 and 6). The two overlapping triangular pulse responses 
h-i(k) and ^(k) shown in FIGS. 5 and 6 then yield the 
desired resulting trapezoidal pulse response hjjk) shown in 
FIG. 4. 

The block diagrams shown in FIGS. 1 and 2 serve only to 15 
illustrate one embodiment of this invention. Elements 1 
through 13 are preferably implemented in the circuitry (as 
hardware). However, implementation in the form of pro- 
gram steps (as software) is also conceivable. 

What is claimed is: 20 

1. A method of estimating a frequency } of a digital input 
signal (x(k)), comprising the steps of: 

determining (2) a phase (4>(k)-arg(x(k)) of the input signal 

(xGO), 

generating (3) a phase difference Ota^iOO) between adja- 25 
cent samples of a phase (<)>(k), <j>(k-l)), 

averaging the phase difference (^^(k)) by using a filter 
(4) which has a trapezoidal pulse response (h^/k)), 
wherein the trapezoidal pulse response (h^/k)) is gen- 3Q 
erated by superimposing a first triangular pulse 
response (h 2 (k)) on a second triangular pulse response 
(h 2 (k)) which is offset in time in comparison with the 
first triangular pulse response (h 1 (k)). 

2. The method according to claim 1, wherein the first 35 
triangular pulse response (h 2 (k)) is generated by folding two 
square-wave pulse responses (h 3 (k), h 4 (k)). 

3. The method according to claim 2, wherein the second 
triangular pulse response (^(k)) is generated by folding the 
first triangular pulse responses (h^k)). 


,684 Bl 

6 

4. The method according to claim 3, wherein the square - 
wave pulse responses (h 3 (k), h 4 (k)) are generated by inte- 
gration (9, 10) and subsequent differentiation (11, 12), where 
a sampling rate (f fll ) is reduced between integration (9) and 
differentiation (11) to one third of the inverse of the period 
(fai^'l/TaJ of the trapezoidal pulse response (h^k)) of 
the filter (4). 

5. The method according to claim 4, wherein the square- 
wave pulse responses (h 3 (k), h 4 (k)) are folded by performing 
the integration (9, 10) and differentiation (11, 12) each in 
two stages. 

6. A device for estimating the frequency (f) of a digital 
input signal (x(k)), comprising: 

a phase determination device (2) for determining a phase 
(<t>(k)=arg(x(k)) of the input signal (x(k)), 

a differentiator (3) for generating a phase difference 
OtaifliM) between adjacent samples of a phase (<(>(k), 
*(k-l)), 

a filter (4) which has a trapezoidal pulse response (h^/k)) 
and averages the phase difference ($ dig Qi)), wherein the 
filter (4) has a first integrator (9), a second integrator 
(10) downstream from the first integrator (9), a sample 
converter (5a) downstream from the second integrator 
(10), a first differentiator (11) downstream from the 
sample converter (5a), a second differentiator (12) 
downstream from the first differentiator (11), and an 
averaging filter (13) downstream from the second dif- 
ferentiator (12). 

7. The device according to claim 6, wherein a sampling 
rate (f al ) is reduced in the sample converter (5a) to one third 
of the inverse of the period (f^'/H/T^ of the trapezoidal 
pulse response (h^k)) of the filter (4). 

8. The device according to claim 6, wherein a transfer 
function H^Z) of the averaging filter (13) amounts to: 

* * * * * 
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